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\^ ' The gluon recombination functions in the twist-4 QCD evolution equations are 

Cp . studied at the leading logarithmic approximation using both the covariant and non- 

covariant methods. We justify that the infrared singularities in the twist-4 coefficient 

functions impede us to simply separate the QCD evolution kernels from the coef- 

o ■ ficient functions at the equivalent particle approximation. In particulary, we point 

JL . out that the gluon recombination functions in the GLR-MQ evolution equation are 

Q I unavailable. The methods avoiding the IR divergences are discussed, which can 

i-Q ■ be used in the derivations of the evolution kernels and coefficient functions at the 



o 



X 



higher twist level. 



1 Introduction 

The QCD evolution equations in the leading order level predict strong rise of parton 
densities when the Bjorken variable xb = Q^f^P-Q decrease toward small values. This 
behavior violates unitarity. Therefore, various models are proposed to modify the twist-2 
evolution kernels (i.e., the parton splitting functions). Gribov, Levin and Ryskin [1] first 
suggest a non-linear evolution equation, in which the evolution kernels (we call them as the 
gluon recombination functions) are constructed by the fan diagrams. Later Mueller and 
Qiu [2] calculate the (real) gluon recombination functions at the double leading logarithmic 
approximation (DLLA) in a covariant perturbation framework. The GLR-MQ equation is 
broadly regarded as a key link from perturbation region to non-perturbation region. This 
equation was generalized to include the contributions from more higher order corrections 
in the Glauber-Mueller formula [3]. Following them, a similar evolution equation (we 
call it as the recombination equation) is derived in a more broad kinematic region (the 
leading logarithmic (Q^) approximation-LLA((5^)) in Ref.[4]. Different from the GLR-MQ 
equation, the gluon recombination functions in this equation are calculated and summed 
in the (old-fashioned) time-ordered perturbation theory (TOPT) [5]. 

In this paper we present the different results of the gluon recombination functions in 
the above mentioned two evolution equations. The cross section (or structure functions) 
of deep inelastic electron-proton scattering in a factorization scheme can be factorized 
into the convolution of coefficient functions with nonperturbative matrix elements, which 
are defined as the parton distributions in the twist-2 level. The QCD evolution kernels 
(for example, splitting and recombination functions) are separated from the coefficient 
functions either in a covariant perturbation theory or in a parton model using TOPT. 



It is well known that these two approaches are completely equivalent. However, we find 
that the infrared (IR) singularities in the twist-4 Feynman diagrams impede us to safely 
extract the evolution kernels before these singularities are cancelled. This result leads 
to the fact that the gluon recombination functions in the GLR-MQ evolution equation 
are unavailable. For illustrating our idea in detail, in Sec. 2 we compare two approaches 
in the derivations of a splinting function in the twist-2 DGLAP equation [6]. Then we 
show two different derivations of the twist-4 recombination functions in Sec. 3. In Sec. 4 
we analyze the IR singularities in the recombination functions, which arise the questions 
in the derivations of the gluon recombination functions in the covariant theory. The 
discussions and conclusions are given in Sec. 5. 



2 Twist-2 splitting functions 

The splitting functions for partons (quarks or gluon) in the twist-2 evolution equations 
are schematically defined as Fig. la. We emphasize that all initial and final partons in 
Fig.l should be on-shell, thus the fraction of the momentum and helicity, which is carried 
by a parton, can be fixed [6]. The physical gauge is convenient in the calculations of 
the evolution kernels. In this work we take the physical axial gauge and let the light-like 
vector n to fix the gauge as n ■ A = 0, A being the gluon field, where 



^^ = ^(1,0^,-1) 



and 



^'^^^(l-Ox,!). (2.1^ 

Therefore, we have 



n ■ n = n ■ n = 0, n ■ n = 1. (2-2) 

Of cause, we can also use other gauges since the result is gauge-invariant. However, 
physical picture of a DIS process is frame- and gauge-dependent. An appropriate choice 
of the coordinate frame and gauge is helpful to extract the evolution kernels using a simple 
partonic picture. 

A natural derivation of the splitting function is given by TOPT in a special infinite 
momentum frame, i.e., the Bjorken frame, where the virtual photon has almost zero energy 
go = (2Mz/ -|- Q'^)/AP and longitudinal momentum q^ = —[2Mv — Q'^)/AP. A Feynman 



propagator is decomposed into the forward and backward propagators in TOPT, where 
the propagating partons stay on-shell, respectively. One can find that the contributions of 
the backward components of the propagators in the cut diagram, for example in Fig. 2a 
are suppressed [4]. Thus, the splitting function can be isolated in the equivalent particle 
approximation [7] as shown in Fig. 2b. To illustrate it, we parameterize the momenta in 
Fig. 2a as 

P= (Po,Pt,P3) = {Pz,0,p,), 



2zpz 



k2 - ((1 - z)p, + ^,. ^^ , , -k^, (1 - z)p,). (2.3) 

2(1 - z)p^ 

The contributions of Fig. 2a to the twist-2 coefficient functions are 
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x§^M{^*q^^*q)^6{z-Xs), (2.4) 

Et Q' 



where we used 



= Tr[-f^kluhTr[-f(,klah (2-5) 

since the momentum k is on-shell. Now we can separately define a splitting function 
Ptopt (Fig- 2b) as 



as dkrp j-,q^q 



1 



SEj^Ek^ [Ej^ + Ek,^ — Ep 
i.e., 



d^k. 



"^(^^^^"'(2^' (2.6) 



P^o'pT - 3 [yZ^) ■ (2-7) 

The same splitting function can be obtained in a covariant framework [8]. The related 
dominant Fynmann diagram for the coefficient function at the LLA{Q'^) is shown in Fig. 
3. For convenience, we take a coUinear (infinite momentum) frame, in which the momenta 
of the initial parton and virtual photon are parameterized as 

P= {Pz,0,p,) =p+n, 

and 

_ g2 

q = q+n + g_n = -xbP+u + n. (2.8) 

2xbP+ 

The contributions of the q —^ qg process in Fig. 3 to the coefficient function are 



C twist— 2 
CVPT 



where 



M{l*q -^ l*q) 



4 1 

= fi-^ < - >coiour -Tr[7^7 ■ ki'jul ■ k'jp'j ■ p'jal ■ k]Ta(3{k2)d'^. (2.10) 

TajS is an axial gauge gluon polarization 



T,Ak)=g,u- ^'\~^^''''r (2.11) 

K ' Tl 

In the calculations we use the Sudakov variables 

k^, = bp^ + cq'^ + k±^, (2.12) 

where two parameters are determined by the on-shell conditions 



b = XB-il-XB)^, (2.13) 

2p ■ q 

and 



c=-^. (2.14) 

2p ■ q 

The result in the leading logarithmic region, where k'^ <^ Q^ are 



ritwist-2 _\^ 2'^s f dk^ 4 1 + Xg I'o 1 c^ 



The splitting function Eq. (2.7) can be directly obtained through dividing the coeffi- 
cient function C^yp^'^ by the contributions of a bare photon-quark vertex (see Fig. 3b), 
which is 

cLe = E ^Mq + p?)M{i*^ - 7*g) 

= Y.el (2.16) 



That is 



^twist—2 /I 1 I ^2 

yq^q ^ ^CVPT _ _Z___ 



ryq^q _ ^ CV FT _ - ^ ' ^ fr, -, j\ 

^CVPT — ~~P^ — o~ T) l^'-^'J 

''bare 



in which we have taken the assumption that 



p+»g_, (2.18) 



thus, one can define 



k+ = xP+. (2.19) 



It is not surprising of the above equivalence of two approaches. In fact, the con- 
tributions of the backward components in two off-shell Feynman propagators in Fig. 3 
are suppressed under the conditions Eq. (2.18) and LLA{Q'^), although straightforward 
computations in covariant theory do not display this fact. For confirming this conclu- 
sion, we recalculate Eq.(2.15) but removing the backward components in two Feynman 
propagators with the momentum k in Fig. 3. We get the same result as Eq. (2.17). 



3 Twist-4 recombination functions 

The recombination functions in the twist-4 evolution equations at LL{Q^)A are de- 
picted in Fig. lb. A simplest way of the derivation of the recombination functions is 
TOPT, which was developed in [4]. In this method the recombination function is isolated 
from the twist-4 coefficient function using the equivalent approximation as shown in Fig. 
4 and it is generally written as 



a^K^OPT dx, ^, - ^g^2 (^^ + ^^)3 \%^P.^kkJ dx, ^^ , [6.1) 

where k"^ = k^ + ky and we use a current (the dashed lines) —^F^^^^Fl^^ to probe the 
gluonic current [8]. 

The momenta of the initial and final partons are parameterized as 



Pi = [xip, 0,0,xip], p2 = [a:2P, 0,0,a;2p], 
k = [xsp + — ^, k±, xsp], k2 = [x4P + 77^^, -k±, X4,p], 

ZXsP 2X4J9 

k^ 

ks = [{X4 - X2)p + —^^, -k_i, {X4 - X2)p], 
2X4P 

k^ 

ki = [{X4 - X2)p + 7-^, -k±, {X4 - X2)p]. (3.2) 

2a;4j9 
For example, the recombination function for GG ^ G at t-channel is 



-TirpQprp — \ — / ■ r^^aXp^l3ri(T^a'X'p'^f3'ri'(T''J'j_ a_|_ 

O color yXi -\- X2) 



rxpik,)Typ,ik4) r,„ _ kikl\ / , _ ^ 



^'""-if^ ^''-^.l (3-3) 



where l,m, r,s are the space indices of a,!], a',!]' of k^ and k4, respectively; GaXp and 
C/3^0- are triple gluon vertex. 
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A set of recombination functions, where the contributions from all channels are summed 
up, are listed in Table I. 

Now we return to discuss the recombination functions in CVPT. The gluon recom- 
bination functions in the GLR-MQ evolution equation are derived using CVPT at the 
double leading logarithmic approximation {DLL A), i.e., at the small x limit [2]. We use 
the same way to recompute the recombination functions but in whole x region. 

The coefficient function containing the gluon recombination function can be written 



as 



d^k 



CS^PT^ = j^,m)m)\MGG-.G?. (3.4) 



In the t-channel. 



\Mgg-.g? = g\-) c"^^c^^"C"'^'^'c^'^'"V^^V'^"'<'rf7' 

O color 

^\p{h)^yp'{ki)^^iy:{ki)T^a{k)Tyai{k)Tr^r^i{k2) 
where the current-gluon vertex v is 



v'^" = k^ki^- g^yk-ki, (3.5) 

and 



/c3 ■ n "*" k^ 

The coefficient function of the contributions of a bare gluonic vertex is 



r.,(A:3) = g., - ^^"t"^^^'" = 4' - 29^nV. (3.6) 



^bare — 9/02" ( -' 



Thus, the recombination function 
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rGG^G 
tdGG^G — ^CVPT /o o\ 

^CVPT — ^G ■ \'^-°) 

bare 

A set of complete recombination functions using the above mentioned CVPT method 
are given in Table II. A surprising result is that the recombination functions with the IR 
singularities in the CVPT method are inconsistent with that in the TOFT method. We 
shall detail it in the following section. 
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4 IR safety in the recombination functions 

At first step, we point out that the IR singularities in Rcvpt originate from the gauge 
term in the gluon propagator T'^^ = g^^ — {k^n^ + k^n^)/k'i ■ n, where k^ ■ n r^ y — x arris 
divergence ii y ^ x. 

We decompose the propagator with the momentum k on the two gluon legs according 
to TOPT in Fig. 5. The on-shell propagating momenta are 



k^ 
k -^ kpizp+n, ^r^^n, k±), (4.1' 

2zp+ 

and 



k^ 
kBi—^-n,zp+n,k±). (4.2) 

2zp+ 

Where the sub-indices F and B refer to the forward and backward propagations, respec- 
tively. The dominant contributions in the LLA^Q"^) are from those terms with the highest 
power of k\ in the numerator of the propagator. Obviously, terms n^ or (!■[ in Eq. (3.6) 
can combine with ks or kp in Eqs. (4.1) and (4.2) through the quark-gluon or triple-gluon 
vertex, respectively. In consequence, the contributions from the forward and backward 
components are coexist in the two propagators with the momentum k in these diagrams. 

It is different from the recombination functions Rcv^pt i ^^e dominant numerator fac- 
tor 7 ■ n in Rqvpt only chooses kp in Eq. (4.2). Therefore, the contributions from 
the backward components are suppressed at the LLA^Q"^) and the equivalent particle 
approximation is available as in Rtopt- 

Such gauge singularities also exist in the derivation of the GLR-MQ equation [2], 
where an unusual ie prescription is used at limit a; ^ 1 and the result is finite after taking 
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the principal value. Obviously, the contributions from the backward components in the 
gluon propagators can not been entirely excluded in this way. Thus, we can conclude 
that the recombination functions in the GLR-MQ equation are really a part of the twist-4 
coefficient functions but not the evolution kernels for parton distribution functions(PDFs) 
since two legs in Fig. 5 are off-shell. 

Of cause, any physical results are IR safe. The IR divergences should be cancelled by 
summing up all related diagrams, including interference and virtual diagrams (see Fig. 
6). After that one can isolate the recombination functions from the coefficient functions 
in the covariant framework. However, at moment there is no available covariant way to 
calculate the virtual and interference diagrams like Fig. 6. In the GLR-MQ equation 
such summations are evaluated using the Abramowsky-Gribove-Kancheli (AGK) cutting 
rule [9] . However, this application of the AGK cutting rule has been argued in Ref . [4] . 
Besides, the AGK cutting rules only change a relative weight in the contributions of the 
real diagrams and they can not cancel any IR divergences. 

Fortunately, the TOPT approach in the Bjorken frame provides an available method 
to derive the recombination functions. In fact, the recombination functions can be safely 
separated from the twist-4 coefficient functions in this way, since the backward compo- 
nents in two parton legs with the momentum k are suppressed. Thus, the equivalent 
parton approximation can be used. Furthermore, the contributions of the gauge singular 
terms disappear due to the absence of these backward components. For justifying this 
conclusion, we use d'^ to replace T'^'^ in Eq. (3.3) and get the same results as Table I. 

In the case that all propagators only keep their forward component in Fig. 6, the 
calculations of the relating interference and virtual diagrams become possible. In fact. 
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a simple relation among the real-, virtual- and interference diagrams is derived in Ref. 
[4]. Thus, we can pick up the contributions of leading recombination from all necessary 
cut diagrams. The resulting evolution kernel for the gluon recombination function has 
following form [4]: 

/ dXidx2f{Xi, X2, Q^)RTOPTi.^l^^2, X b) 

J2{x\+X2)>xb 

-2 / dXidX2f{Xi,X2,Q^)RT$PT{^l^^2,XB), (4.3) 

J{x\+X2)>XB 

where f{xi,X2,Q^) is the two gluon correlation function. The first and second terms in 
Eq. (4.3) give the antishadowing and shadowing effects in the evolutions of the parton 
densities. 

Now let us look back to the CVPT method for the derivation of the recombination 
functions. According to the above mentioned discussions, we use (i_L to instead Txp in Eq. 
(3.6). As expected, one can get the same results as Table I. From the above mentioned 
discussions we can conclude that the forward propagators dominate the the recombination 
functions. 
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5 Discussions and summary 

An interesting question is why the IR singularities of the gluon propagator in the twist- 
2 coefficient function (Fig. 3a) does not break the equivalent particle approximation. We 
note that this propagator through the dashed line is on-shell. Thus, Eq. (2.11) can be 
rewritten as 

f"^ = df - ^n'^n" (5.1), 

Now the dominant contributions to the LL{Q'^)A are only from the forward components 
and the equivalent particle approximation is applicable. 

Our discussions for the recombination functions can also be used in the investigation of 
the twist-4 coefficient functions [11]. One can get those functions up to dk\/k\ accuracy 
by the product of the recombination functions in Table I and the bare vertex (2.16) or 
(3.7). 

In summary, we show that the QCD evolution kernels are frame- and gauge-independent, 
however, the separation of these kernels from the coefficient functions depend on the frame 
and gauge. The IR singularities in the recombination functions inhibit us to safely isolate 
the recombination functions. This result leads to the fact that the gluon recombination 
functions in the GLR-MQ evolution equation are unavailable. The methods avoiding the 
IR divergences are discussed by using the TOPT. 
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Table. 1 



oGG-^qq 

rirpQpj, 


1 {2y-xf[l8y^-21yx+Ux^) 
48 ?y5 


rtqq^qq 
rtrpoPT 


2i2y-xr 
9 y^ 


T3qG^qG 
^TOPT 


1 {2y-x)[l6y^+20yx+25x^) 

72 v** 


oqq^GG 
JXtoPT 


8 (2y^x)[-10y'^+5yx+2x^) 
27 y'-^x 


jdGG^GG 


9 {2j/-a;)(72j/4~48j/'''a;+140j/2a;2- 


-116j/a;^+29a;4] 


64 ?y5^ 1 





Table. 2 




TiGG^qq 
^CVPT 


1 (2y-x)2(l8j/2-21yx+14x2) 
48 v^ 


jjqq^qq 
^CVPT 


2 (2y-x)^(l3j/2+6ya;+3a;2] 
27 vHy-x? 


T^qG^qG 
^CVPT 


I (2 ?;-x) (85 y-'-gS y^x+ioe y^^^a.gg ^,3,3+3]^ ^,4 j 




36 yHv-x)' 


jjqq^GG 
^CVPT 


g (2j/-a:)(-10j/2+5j/3:+2a:2) 
27 y^x 


pGG^GG 


9 (2 jz-x) ( 144 / -360 j/5a;+54i j^4^2 _^^q y'-^x^+22l y'^x'^- 


-70j/a;5+12a;6] 


64 j/5(j/-x)''a; 
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Figure Captions 

Fig. 1 The schematic descriptions of (a) sphtting and (b) recombination functions. The 
sohd hnes are gluon or quarks and the grey ovals imply all possible QCD channels. 

Fig. 2 Graphical illustration of the isolation of a splitting function in the TOPT cut 
diagram, k is on-shell momentum and dashed lines refer the time order. 

Fig. 3 (a) A dominate covariant Feynman diagram for the twist-2 coefficient function at 
the LLA; (b) the bare virtual photon-quark vertex. 

Fig. 4 Graphical illustration of the isolation of a recombination function in the TOPT 
cut diagram.The dashed line is an effective "current" probing the gluonic matrix. 

Fig. 5 A covariant Feynman diagram containing the gluon recombination function. 
Fig. 6 The virtual and interferant diagrams relating to the gluon recombination function. 
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